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We analyze the decay modes of D/Ds — >■ PP, PV on the basis of a hybrid method with the 
generalized factorization approach for emission diagrams and the pole dominance model for the 
annihilation type contributions. Our results of PV final states are better than the previous 
method, while the results of PP final states are comparable with previous diagrammatic 
approach. 



1 Introduction 

The CLEO-c and the two B factories already give more measurements of charmed meson decays 
than ever. The BESIII and super B factories are going to give even much more data soon. 
Therefore, it is a good chance to further study the nonleptonic two-body D decays. However, 
it is theoreticahy unsatisfied since some model calculations, such as QCD sum rules or Lattice 
QCD, are ultimate tools but formidable tasks. In B physics, there are QCD-inspired approaches 
for hadronic decays, such as the perturbative QCD approach (pQCD)IJJ the QCD factorization 
approach (QCDF)E] and the soft-collinear effective theory (SCET)P But it doesn't make much 
sense to apply these approaches to charm decays, since the mass of charm quark, of order 1.5 
GeV, is neither heavy enough for a sensible l/rric expansion, nor light enough for the application 
of chiral perturbation theory. 

After decades of studies, the factorization approach is still an effective way to investigate 
the hadronic D decaysSl. However, the naive factorization encounters well-known problems: the 
Wilson coefficients are renormalization scale and 75-scheme dependent, and the color-suppressed 
processes are not well predicted due to the smallness of 02. The generalized factorization ap- 
proaches were proposed to solve these problems, considering the significant nonfactorizable con- 
tributions in the effective Wilson coefficients^. Besides, in the naive or generalized factorization 
approaches, there are no strong phases between different amplitudes, which are demonstrated 
to be existing by experiments. 

On the other hand, the hadronic picture description of non-leptonic weak decays has a longer 
history, because of their non-perturbative feature. Based on the idea of the vector dominance, 
which is discussed on strange particle decays^ the pole-dominance model of two-body hadronic 
decays was proposed^ This mo del h as already been applied to the two-body nonleptonic decays 
of charmed and bottom mesons^^^. 

In this work, the two-body hadronic charm decays are analyzed based on a hybrid method 
with the generalized factorization approach for emission diagrams and the pole dominance model 
for the annihilation type contributions^. 



2 The hybrid method 



In charm decays, we start with the weak effective Hamihonian for the AC = 1 transition 

Gf 

-Heff = -^VcKuiCiOi + C2O2) + h.c, (1) 

with the current-current operators 

01 = nQ7^(l - 75)^2/3 • 93/37^(1 - 75)Ca, 

02 = Ua7^(l - 75)^20 • 93/37''(l - 75)C/3. (2) 
In the generahzed factorization method, the amplitudes are separated into two parts 

{MiM2\neff\D) = -^FcxMai,2(Mi|gi7^(l - 75)9210) (M2|g37''(l - l5)c\D), (3) 

where ai and 02 correspond to the color-favored tree diagram (T) and the color-suppressed 
diagram (C) respectively. To include the significant non-factorizable contributions, we take 01^2 
as scale- and process-independent parameters fitted from experimental data. Besides, a large 
relative strong phase between ai and 02 is demonstrated by experiments. Theoretically, the 
existence of large phase is reasonable for the importance of inelastic final state interactions in 
the charmed meson decays, with on-shell intermediate states. Therefore, we take 

ai = |ai|, a2 = |a2|e*'', (4) 

where ai is set to be real for convenience. 

On the other hand, annihilation type contributions are neglected in the factorization ap- 
proach. However, the weak annihilation (I^-exchange and H^-annihilation) contributions are 
sizable, of order l/rric, and have to be considered. It is also demonstrated to be important by 
the difference of life time between and D~^. The pole-dominance model is a useful tool to 
calculate the considerable resonant effects of annihilation diagrams. For simplicity, only the 
lowest-lying pole is considered in the single-pole model. Taking — t- PP, PV as example, 
the annihilation type diagram in the pole model is shown in Fig(Ha). goes into the inter- 
mediate state M via the effective weak Hamiltonian in Eq.([T]), shown by the quark line in the 
FiglU^b), and then decays into PP{PV) through strong interactions. Angular momentum should 
be conserved at the weak vertex, and all conservation laws be preserved at the strong vertex. 
Therefore, the intermediate particles are scalar mesons for PP modes and pseudoscalar mesons 
for PV modes. In decays, they are IF-exchange diagrams, but T^-annihilation amplitudes 
in the decay modes. 
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Figure 1: Annihilation diagram in the pole-dominance model 
The weak matrix elements are evaluated in the vacuum insertion approximatioiP, 
{M\'H\D) = -^VcKMaA,Emqil^.{l-lM^)m^l^{^-lMD) 

= ^VcKMaA,EfMfDml), (5) 



where the effective coefficients a a and qe correspond to PF-annihilation and Vl^-exchange amph- 
tudes respectively. Strong phases relative to the emission diagrams are also considered in these 
coefficients. 

For the PV modes, the effective strong coupling constants are defined through the La- 
gran gi an 

£ypp = igvPpV^'{Pid^P2 - P2d^Pi), (6) 

where gvpp is dimensionless and obtained from experiments. By inserting the propagator of 
the intermediate state M, the annihilation amplitudes are 

Gp 1 

(PVinefflD) = —7=VcKMaA,EfMfDm\,— ^gvPM'^ie* -Pd)- (7) 

y/ 2 mj-,-mj^,j 

As for the PP modes, the intermediate mesons are scalar particles. The effective strong coupling 
constants are described by 

J-SPP = -gsppmsSPP. (8) 

However, the decay constants of scalar mesons are very small, which is shown in the following 
relation 

fs rns 

where fs is the vector decay constant used in the pole model, fs is the scale-independent scalar 
decay constant, mi^2 are the running current quark mass, and ms is the mass of scalar meson. 
Therefore, the scalar pole contribution is very small, resulting in little resonant effect of annihi- 
lation type contributions in the PP modes. On the contrary, large annihilation contributions are 
given in the PV modes by relative large decay constants of intermediate pseudoscalar mesons. 



(9) 



3 Numerical results and discussions 

In this method, only the effective Wilson coefficients with relative strong phases are free param- 
eters, which are chosen to obtain the suitable results consistent with experimental data. For 
PP modes, 



Ol = 0.94 ±0.10, 02 = (0.65±0.10)e*(^^2^^°)° 

aA = (0.20±0.10)e^(300±io)°^ a^, = (1.7 ± 0.1)e*(90±io)° _ (^q) 



For PV modes, 



ai 



.1 = 1.32 ±0.10, ar = (0.75±0.10)e»i°)°, 

a^/ = (0.12±0.10)e^(345±ior, ^ = (0.62 ± 0.10)e^(238±io)° _ (j^) 

All the predictions of the 100 channels are shown in the tables of refP. The prediction of 
branching ratio of the pure annihilation process Df — )• vr'^vr'^ vanishes in the pole model within 
the isospin symmetry. It is also zero in the diagrammatic approach in the flavor SU(3) symmetry. 
Simply, two pions can form an isospin 0,1,2 state, but is ruled out because of charged final 
states, and isospin-2 is forbidden for the leading order AC = 1 weak decay. The only left s-wave 
isospin-1 sate is forbidden by Bose-Einstein statics. In the pole model language, G parity is 
violated in the isospin-1 case. Therefore, no annihilation amplitude contributes to this mode. 

The theoretical analysis in the r] — rj' sector is kind of complicated. The predictions with r]' 
in the final state are always smaller in this hybrid method than those case of ij due to the smaller 
phase space. However, it is opposite by experiments in some modes, such as — )■ 7r+r/(77'), 



— )■ K^ri{rj'). This may be the effects of SU(3) flavor symmetry breaking for riq and rjs, the 
error mixing ang le between , and inelastic final state interaction, or the two gFuon anomaly 
mostly associated to the 77', etc.. The mode of — p'^r]{r]') is similar with the above two 
cases, the opposite ratio of r] over rj' between theoretical prediction and the data. But this is a 
puzzle b y ex periment measurement, which is taken more than ten years ag As is questioned 
by PDG^^, this branching ratio of (12.5 it 2.2)% considerably exceeds the recent inclusive ij' 
fraction of (11.7 ± 1.8)%. 

Recently, model independent diagrammatic approach is used to analyze the charm decays^^^l. 
All two-body hadronic decays of D mesons can be expressed in terms of some distinct topological 
diagrams within the SU(3) flavor symmetry, by extracting the topological amplitudes from the 
data^^. Since the recent measurements of — )• Tr^p^^^ and — )• tt^w^ISI g^ye a strong 
constraint on the TV— annihilation amplitudes, one cannot find a nice fit for Ap and Ay in the 
diagrammatic approach to the data with Df — t- K*^K^,K^K*^ simultaneously. Compared to 
the calculations in the model-independent diagrammatic approach ^1^, our hybrid method gives 
more predictions for the PV modes in which the predictions are consistent with th e exp erimental 
data. It is questioned that the measurement of Br[Df — t- K'^K*^) = (5.4 it 1.2)%^^^ which was 
taken two decades ago, was overestimated. Since |Cy| < \Cp\ and Ay ~ Ap as a consequence of 
very smah rate of 7r+/9°, it is expected that Br{D+ -f K^K*+) < Br{D^ K*^K+) = 

(3.90 lb 0.23)%. Our result in the hybrid method also agrees with this argument. 

As an application of the diagrammatic approach, the mixing parameters x = (mi — m2)/V 
and y = {Ti — r2)/r in the — mixing are evaluated from the long distance contributions 
of the PP and VP mo destlSl. The global fit and predictions in the diagrammatic approach are 
done in the SU(3) symmetry limit. However, as we know, the nonzero values of x and y come 
from the SU(3) breaking effect. Part of the fiavor SU(3) breaking effects are considered in the 
factorization method and in the pole model. Therefore, our hybrid method takes its advantage 
in the analysis of — mixing. 
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